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ABSTRACT 



The transit of Venus in 2004 offered the rare possibility to remotely sense a well-known planetary atmosphere using ground-based 
observations for absorption spectroscopy. Transmission spectra of Venus' atmosphere were obtained in the near infrared using the 
Vacuum Tower Telescope (VTT) in Tenerife. Since the instrument was designed to measure the very bright photosphere of the Sun, 
extracting Venus' atmosphere was challenging. CO2 absorption lines could be identified in the upper Venus atmosphere. Moreover, 
the relative abundance of the three most abundant CO2 isotopologues could be determined. The observations resolved Venus' limb, 
showing Doppler- shifted absorption lines that are probably caused by high-altitude winds. 

This paper illustrates the ability of ground-based measurements to examine atmospheric constituents of a terrestrial planet atmosphere 
which might be applied in future to terrestrial extrasolar planets. 
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1. Introduction 

Transits of Venus, during which Venus crosses the line of sight 
between Earth and the Sun, come in pairs separated by eight 
years over a total period of ~121.5 years. The last transit oc- 
curred in 2004 whereas the next one will be observable in 2012. 
Apart from the inner planets of the solar system, extrasolar plan- 
ets also have a certain possibility to transit in front of their cen- 
tral star, offering the feasibility of investigating their atmosphere 
during their primary transit. Up to now, 140 transiting extrasolar 
planets are known out of the more than 500 exoplanets that have 
been found in totaQ. Out of this, seven transiting pl anets might 



be cla ssified as terrestrial: Co RoT-7 b (4.80MEart h; Leger et all 
I2009h . GJ 1214 b (6.36MB,„t.. lCharbonneau et al.ll2009l) as well 
as the recently detected plan ets Kepler- 11 b,d,e,f (4.30, 6.10, 
8.40, 2.30MF^rth. respectivelv.[Lissauer et al.ll201 lb and Kepler- 
10 b (4.54Mp.,.h. rBatalha et al.ll201 ih . 

During a planetary transit, the fight emitted by the central 
star traverses the planetary atmosphere and is attenuated before 
it reaches the observer Analysis of the extinction allows to char- 
acterize the planet, for example the chemical composition of the 
atmosphere. The apparent planetary radius at which the tangen- 
tial optical depth is below unity varies at different wavelengths 
and can hence be used to infer the height of the atmosphere (re- 



' See J. Schneider's Extrasolar Planet Encyclopedia at 
|http://exoplanet.eU) (status: June 201 1) 



lated to pressure and atmospheric structure) and the existenc e of 
cloud layers (see e.g. lSeager & Sasselovll2000tlBrownll200lh . 

In the past, transits of Mercury offered the opportunity to de- 
tect absorption features within the atmosphere by ground-based 
spectroscopic observations. For example, during the Mercury 
transit in 2003, neutral sodi um in its exosphere ha s been de- 
tected using this technique dSchleicher et al.l l2004l) . Also the 
Earth can be investigated as a transiting planet during lunar 
eclipse observations, showing absorption features of the major 
atmospheric constituents like CH4, H2O, O2, O3 a nd CO? in 
reflected light from the Sun by the Moon (see e.g. iPalle et aTl 
I2OO9I : lyidal-Madjar et al.ll2010h . Furthermore Earthshine obser- 
vations, where the reflected sunlight from Earth is observed on 
the Moon's night side, show absorption signatures of O2, O3 
and H2O as well a s a s ignature of Earth's vegetation red edge 
dArnoId et aI.ll2002l:lAmold 2008i) . 

This technique can of course also be applied to transit- 
ing exoplanets: Molecular absorption bands in the infrared 
have been dete cted in the atmosphere of planets around 
brigh t stars (e.g. lKnutsonll2007l: iTinetti et alJl2007l: ISwain et all 
l2008h . Furthermore, the sodium resonance doublet at 598.3 nm 
in the atmosph ere of the transiting ex trasolar giant planet 
HD 209458 b dCharbonneau et al.1 l2002h during four plane- 
tary transits has been detecte d, as well as neutral hydrogen 
(HI) ( Vidal-Madi ar et al.l|2003h . oxy gen (OI) and c arboii (CII) 
dWlal-Madiar et al.ll2004l) . Recendv. ISnellen et alJ (1201 Oh used 
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high-resolution ground-based spectroscopy to detect CO hues in 
the atmosphere of HD 209458 b at about 2 300 nm. These Unes 
showed a blue shift of about 2 km s ' with respect to the veloc- 
ity of the host star, perhaps indicating a velocity flow from the 
day to the night side, also known as subsolar-to-antisolar (SSAS) 
flow. 

However, performing transmission spectroscopy of terres- 
trial exoplanets and detecting the faint atmospheric absorption 
lines is a difficult task. Up to now, only one terrestrial ex- 
oplanet has been cha racterized spectro scopically using trans- 
mission spectroscopy: iBean et al.l (1201 Oh obtained transmission 
spectra from the terrestrial exoplanet GJ 1214b in the range from 
780 to 1 000 nm, showing a lack of spectral features which in- 
dicates either a dense, water vapor dominated atmosphere or a 
hydrogen dominated atmosphere with a cloud layer which is 
opaque i n the obser v ed wa velen gth range. Addition al observa- 
tions by ICroll et al.l (1201 ll) and (Pesert et all (1201 lb rule out a 
fully water dominated atmosphere. Currently a solar composi- 
tion atmosphere witho ut methane and Rayleigh scatterin g by 
clouds best fits the data dMiller-Ricci Kempton et alj|201 Ih . 

Major modeling efforts for synthetic spectra of hy- 
pothetical terre s trial e xoplanets were perform e d by 
iDes Marais eTail (|2002|): ISegura et al.l (I2003L l2005l) : iTinettil 



(120061): lEhrenreichetal.1 ()2006h 



- . iKaltenegger et al.l (l2007l) . 

iKaltenegger & Traubl (120091) and iRauer et all (1201 ll) . in which 



the terrestrial planets of the solar system (i.e. Venus, Mars and 
Earth) are investigated as if they were distant exoplanets. These 
studies also include different parameter studies, in which e.g. the 
central star was varied , the orbital distance or the atmospheric 
composition. Note that lEhrenreich et all (|2006|) has studied how 
Venus would look like if seen as a transiting exoplanet, when 
orbiting different main-sequence stars. 

This paper investigates the feasibility of the transmission 
spectroscopy technique for analyzing the atmosphere of a terres- 
trial (exo-) planet by using ground-based high-resolution mea- 
surements during a Venus transit. A cloud-free radiative trans- 
fer model was used to compare these measurements with syn- 
thetic transmission spectra of Venus' upper atmosphere above 
the cloud top, assuming a pure CO2 atmosphere. 

In Sect. |2] the observations and the data reduction are de- 
scribed and the extracted limb spectra are discussed. The radia- 
tive transfer model and the calculations are described in Sect.|3] 
In Sect, m a parameter variation is performed, before the final fit 
to measured data is performed in Sect. |5] Finally the results are 
summarized in Sect.|6]and implications for extrasolar planets are 
discussed in Sect.|7] 

2. Observations and Reduction 

The transit of Venus on June 8th, 2004 was observed at the 
Vacuum Tower Telescope (VTT) of the Kiepenheuer Institut fiir 
Sonnenphysik (KIS) in Tenerife. The observations of the tran- 
sit lasted from 06:30h (UT), when the Sun rose with Venus al- 
ready in front of the disk, until ll:26h (UT) at fourth contact. 
The VTT is a 70 cm diameter telescope with a focal length of 
46 m, equipped with an Echelle spectrograph. The spectrograph 
slit corresponds to a width of 0.35" and a length of 92" on the 
sky. In the focal plane of the spectrogr aph, the Tenerife Infrared 
Polarimeter (TIR ICollados et akll 19991) detector was installed. It 
consists of four 128x128 pixel CCDs and was used in its spec- 
trograph mode (without polarimeter). The spectral resolution of 
the entire system was R - A/AA = 200 000. The image scale in 
the focal plane was 0.35" px'. Thus one pixel corresponds to 
about 71 km at the position of Venus. 



Three wavelength regions (named by the observers) with 
strong absorption lines of CO2 were observed with a sampling 
of0.031Apx-': 

- 1 596.49nm - 1 597.29 nm ("Favorite") (1 278 images) 

- 1 597.49 nm - 1 598.27 nm ("Favorite-i-") (4779 images) 

- 1612.39nm- 1 613.19nm ("Grabbag") (4 103 images) 

One single image is an integration of 20 exposures with an 
exposure time of 50 ms to provide 1 s integration time. Under the 
same conditions, 207 flatfield images on the unobscured solar 
disk and 127 darks were taken. 

The telescope tracked the Sun and drift scans over the Venus 
disk were performed, i.e. Venus moved across the slit that was 
fixed with respect to the solar image. For the observations in the 
"Favorite" wavelength range, the slit was oriented at an angle 
of 315.5° (measured from Venus' north pole, anti-clockwise). 
During the observations in the "Grabbag", as well as in the 
"Favorite+" wavelength range, the slit was rotated to an angle of 
45.5° (measured from the Venusian north, anti-clockwise) and 
the observati on was further a ccomplished within this position. 

Note that lAlonsol (|2006"} examined the data with respect to 
latitudinal variations of absorption line positions, observing dif- 
ferent depths for the spectral features in the "Fayorite+" region 
at higher and lower latitudes, but without further interpretation 
of the data. In this paper, however, an average over the different 
latitudes is used and the absorption of the three most abundant 
isotopologues of CO2 is thus investigated with a higher signal- 
to-noise ratio (SNR). 

Note that the "Grabbag" region finally turned out not to be 
usable due to a strong noise pattern and a strong gain jump be- 
tween the right and left CCDs, which could not be corrected. 
Therefore transmission spectra of this region are omitted in the 
following. 

An example of a raw image (Is integration time) is shown 
in Fig.[T^. It shows a ring-like fringe pattern, a central horizon- 
tal line separating the two upper 128x128 pixel CCDs from the 
lower ones, and a region with a noticeable number of bad pix- 
els in the upper left corner. In the figure the y-axis is the posi- 
tion, whereas the x-axis is the wavelength. Slightly dark vertical 
bands show the solar and terrestrial absorption lines. The disk of 
Venus can be seen as a thick horizontal band in the center of the 
image. 

A standard image reduction procedure has been used to cali- 
brate the images, including dark current removal and flatfielding. 
Note that the final flatfield master image did not contain any ab- 
sorption lines of the Sun or Earth, but only the fringe pattern. 
Thus in the final reduced images (see Fig.[T]3), only the interfer- 
ence fringes are removed whereas solar and telluric absorption 
lines are now clearly visible as well as the separation between 
the CCDs. 

The next step was to extract the spectra of Venus' atmo- 
sphere. First, the position of both Venus limbs were identified 
on the reduced images and taken as reference points. The limb 
positions were determined by fitting the peaks of the first deriva- 
tive of the spatial intensity with a Gaussian. 

The mean diameter (taking into account the pixel scale of the 
CCD) of the projected Venus disk in the "Favorite" region was 
7 980.018±2 164.985 km, whereas in the "Favorite+" region it 
was 4 683. 8790+97 1.466 km. Thus the slit was never positioned 
over the center of the Venus disk, but more to the limb, and 
hence the lower limit on the projected atmospheric altitude cov- 
ered by one pixel is 46. 313+12. 565 km in the "Favorite" region 
and 27. 184+5. 638km in the "Favorite+" region, assuming that 
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Fig. 1. a) Raw image from the "Favorite+" region: Clearly visi- 
ble is the Venus disk as the horizontal dark strip in the center of 
the image, the fringe pattern (sloping dark lines), the bad pixel 
area in the upper left corner and the gain jump between the up- 
per and lower CCD chips in the center of the frame. Strong solar 
and weak telluric absorption lines appear as vertical dark stripes, 
b) Reduced image, after flatfielding. Clearly visible are the so- 
lar and telluric absorption lines and the gain jump between the 
CCDs. 



the radius of the visible disk of Venus is 6 1 16.80km (hence in- 
cluding the opaque cloud deck at 65 km altitude). Therefore the 
upper atmosphere can be resolved within a few pixels. 

After that, a spectral and spatial normalization was per- 
formed: First, each image row was divided by a spectrum of the 
solar disk, which was obtained by binning two image rows 1 1 
pixels offset from the Venus limb. This step removes the telluric 
absorption lines but only attenuates the solar absorption lines. 
Note that in order to distinguish between telluric absorption lines 
and lines forming in the atmosphere of an exoplanet, ground- 
based observations of transiting exoplanets may take advantage 



of the Doppler shift due to the prop er motions of the star, th e 
exoplanet, the Sun and Earth (see e.g. lVidal-Madiaret al.ll2010l) . 

Next, the spatial normalization was performed by averaging 
image columns in a clean part of the CCD and dividing the image 
by this. The resulting image in Fig. |2] now clearly shows weak 
absorption lines in the transition region between the Venusian 
and the solar disk, which do not correlate with the solar absorp- 
tion lines. Comparison with the HITRAN2004 ( High-Resolution 
Trans mission Molecular Absorption Database, iRothman et al.l 
l2005h Une database identifies the atmospheric absorption lines 
as those of CO2. The flatfielding did not work properly at the 
position of the Venus disk, as the fringe pattern is still visible. 

The resulting residuals of solar absorption lines visible in 
Fig. |2] arise from Doppler shifts of convecting granules in the 
photosphere of the Sun and the so called "5 minute" oscillation. 
The Doppler shift of these residuals is ~0.0143 nm, which corre- 
sponds to a velocity of 1 .7 km s ' and a size of ~1 300 km. These 
values are consistent with values measured by ICarroU & Ostlid 
(Il996h . 



Solar disk 



Atmospheric layer 



Absorption lines 



Residuals of solar granulation 



Fig. 2. Flatfielded and normalized image from the "Favorite+" 
region. In the transition region between the solar surface and the 
disk of Venus a thin layer with absorption features is visible, 
which represents the atmospheric layers of Venus. 

Finally, the spectra of Venus' atmosphere at the limb position 
were extracted and all spectra of a particular wavelength region 
were averaged, in order to improve the SNR. Thus the informa- 
tion about the altitude and the corresponding Doppler-shifts by 
winds are lost. Note that the afore mentioned normalization is 
not perfect, primarily due to varying gain jumps between the 
CCDs and at wavelengths with strong solar absorption lines. 
This results in a distorted baseline in the averaged spectrum at 
each limb. The baseline in every image is therefore corrected 
by applying a ninth degree polynomial fit. It was checked that 
the relative absorption depths were not affected by this correc- 
tion. Furthermore, on some images the atmospheric layers are 
disturbed by the fringe pattern within the Venus disk, inhibiting 
any detection of absorption lines. These images were removed 
manually before averaging all spectra. 

In total 395 images from the "Favorite" and 1 107 images 
from the "Favorite+" region have been averaged. Using line po- 
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sitions from HITRAN as a reference, an accurate wavelength 
calibration was performed. See Figs.[3]a) and b) for the averaged 
spectra from the "Favorite+" and "Favorite" wavelength range, 
respectively. Note that the Doppler shift of absorption lines due 
to the high wind velocities in the upper atmosphere of Venus was 
neglected, although superrotation plays a major role in the upper 
atmosphere dynamics. 

However, still a Doppler shift between adjacent pixel rows 
can be measured: Note that the Venus absorption lines in Fig.|2] 
are tilted on both sides of the limb, with a blue shift of the lines 
with increasing altitude on the eastern Venus limb (lower part of 
the Venus disk in Fig. |2]i, and a red shift of the absorption lines 
on the western limb. The Doppler shifted lines yield a velocity 
gradient of about 7ms' per km altitude, which remains nearly 
constant with increasing altitude on the western limb, while it 
decreases slightly on the eastern limb. 

This might be caused by a combined effect of the superrota- 
tion of Venus' mesosphere an d the SS AS flow in the upper atmo- 
sp here ([Bougher et al. 19970 . SS AS wind velo cities measured 
bv [Goldstein et alj ( ll991h and lBetz et al.1 d 1977b at an altitude of 
100 km are in th e range from 105 to 1 40 ms~', whereas recent 
measurements bv lLellouch et al.l (l2008h from the Venus Express 
mission show wind speeds in the range from 30 to 50ms ' 
at -93 km altitude strongly increasing up to 90 to 120 ms"' 
at ~ 102 km altitude. The velocity gradient measured from the 
Venus trans it spectra is thus in ag reement with the gradient de- 
termined bv lLellouch et al.l (l2008l) . A detailed study of this effect 
is however beyond the scope of this paper and will be the subject 
of a following paper, perhaps incorporating observations during 
the next Venus transit in June, 2012. 



3. Radiative transfer calculations 

The radiative transfer model SQuIRRL (Schwarzschild 

Quadrature InfraRed Radiation Line-by-line, 

ISchreier & Schimpfl 1200 lb has been used to fit the mea- 
sured spectra. SQuIRRL is a cloud free, line-by-line radiative 
transfer model which uses the molecular absorption line 
database HITRAN for the calculation of absorption cross 
sections. SQuIRRL calculates emission and absorption spectra 
for arbitrary atmospheric paths, assuming local thermodynamic 
equilibrium (LTE). Furthermore a convolution with appropriate 
spatial and spectral instrumental response functions can be 
performed. 

The model atmosphere extends from to 180 km altitude, 
divided into 30 layers. CO2 is treated as the only atmospheric 
absorber with an isoprofile of 96.5% abunda nce up to 10 0km 
altitude. Above, the CO2 profile from iHedin et al.l (Il983h was 
used. For the transmission calculations, below 100 km altitude 



the day side temperature profile fror n ISeiff et al.l (11985b was 
used, whereas above the profile from iHedin et aP d 1983b was 
adopted, which is in agreement with recent solar occultation 
measurements o f Venus' mesosphere from the Venus Express 
mission (see e.g. lMahieux et al.ll2010b . 

Transmission spectra for a set of different tangential heights 
above 100 km altitude are calculated for each wavelength region, 
using a Gaussian beam of 24 km and 15 km half-width half- 
maximum (HWHM) for the "Favorite" and "Favorite-H" range, 
respectively, according to the average projected atmospheric 
height in the respective wavelength region. Note that although 
one pixel corresponds to about 71 km at the position of Venus as 
mentioned in Sect. |2l the slit was not positioned perpendicular 
to the surface, but under a certain angle since the spectra were 




a) 



1597.8 1598.0 1598.2 

Wavelength [nm] 




b) 



1596,8 1597,0 

Wavelength [nm] 



Fig. 3. Averaged spectrum (solid line) of both Venus limbs in 
the a) "Favorite+" and b) "Favorite" region together with the 
best fitting model spectrum (dash-dotted line, see text for de- 
tails). Absorption line positions of the three isotopologues are 
indicated by different symbols. 

The dashed line shows a scaled spectrum of the Sun in arbitrary 
units, showing partial overlap of solar absorption fines with ab- 
sorption lines from Venus' atmosphere. 



taken close to the limb of the Venus' disk. Thus the projected at- 
mospheric height is used in order to account for the field of view. 
The resulting spectra are finally convolved with a Gaussian of 
6.15 X 10"^ nm HWHM to account for the spectral resolution of 
the instrument. 

Note that in the radiative transfer model no clouds or aerosol 
opacity is included. Nevertheless, the measured spectra can be 
fitted quite well since the best fitting tangential heights of the 
extracted spectra are well above an altitude of 100 km, as will 
be shown later Above 100 km, clouds can be neglected to a 
good approximation since extinction by H2SO4 cloud particles 
is expected to extend only up to abou t 80 km altitude (see e.g . 
lEspositoetalll 19831: iRoos et al.|[T993l or lGrinspoon et alj|1993b 
Venus Express SPICAV/SOIR observations by IWilquet et al.l 
( l2009b have shown that the extinction coefficient of the atmo- 
spheric haze measured at 1 553.7 nm decreases by about two or- 
ders of magnitude to about 2 x 10"'* km"' at 90 km, which is by 
far lower than the CO2 absorption coefficient and thus can be 
neglected for the altitudes probed in this paper. 

Furthermore, the assumption of LTE for the radiative trans- 
fer calculati ons may not hold in the altitude range con s idered 
in this paper. lRoldanetalJ ( l2000b : lL6pez-Valverde et alj ( l2007b 
have found a departure from LTE conditions above 90 km al- 
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titude (depending on the CO2 band), producing partly strong 
emission lines. Nevertheless, non-LTE conditions in the wings 
of the 1 600 nm absorption band, as considered in this paper, can 
be neglected due to several reasons: 

- Under non-LTE conditions, the excitation temperature of the 
CO2 (300 11, 30012, 30013) levels (which correspond to the 
observed absorption lines) by solar excitation is about 400 K 
at altitudes above about 120 km in the wavelength range con- 
sidered here (M. Lopez-Puertas, priv. comm.). This temper- 
ature is higher than the atmospheric temperature of about 
300 K under LTE conditions (see Fig.|4]i. However, the solar 
flux with a blackbody temperature of about 5 800 K is much 
higher than the flux emitted in the Venus atmosphere. 

- iLopez-Valverde et al.l (l2007h found only a very weak atmo- 
spheric emission (which is below the detection limit of the 
VIRTIS instrument on Venus Express) for the 1 600 nm CO2 
band at the limb for low solar zenith angles (i.e. for the trans- 
mission geometry considered in this paper) which is even 
lower in the band wings where the Venus transit observations 
have been performed. 

- Potential non-LTE effects caused by the non-LTE popula- 
tions of the lower vibrational level of hot bands can also be 
neglected, since the observed lines correspond to a funda- 
mental transition from the ground vibrational level. 

To summarize, the contribution of the non-LTE atmospheric 
emission is in any case much smaller than the Sun's radiation 
and thus strong lines seen in emission are not expected in the 
considered wavelength range. Also an influence on the absorp- 
tion can be neglected. Note that observations of the next Venus 
transit in 2012 may provide a much better data quality and may 
be achieved at wavelengths where non-LTE effects can not be ne- 
glected. In that case a more detailed radiative transfer code needs 
to be used, which includes non-LTE effects as well as cloud and 
aerosol extinction. 

Note that for a transiting exoplanet the temperature pro- 
file is not available. The day side temperature profile might 
be inferred from atmospheric modeling and secondary eclipse 
observations by investigating the emission spectrum. However, 
for transmission observations the terminator temperature profile 
cannot be inferred easily, since it depends on numerous factors, 
like planetary rotation, atmospheric composition and dynamics, 
distance to star, etc.. Nevertheless, trans mission obs e rvatio ns of 
the Hot- Jupiter planet HD 209458 b by ISing et al.l (I2OO8I) and 
[Vidal-Madiar et al. (2011]) were able to infer a temperature pro- 
file at the terminator by investigating sodium absorption line pro- 
files. Close-in exoplanets which may have a thick atmosphere 
and rotate slowly will be more likely to have a strong SSAS flow 
in order to transport warm air from their hot day side to the cold 
night side, thus a day side temperature profile might be appli- 
cable for the terminator region. Furthermore these atmospheres 
will have a higher atmospheric s cale height on the day side (see 
for example burrows et al.l2010i showing this effect on a transit- 
ing giant exoplanet), which could in principle be detected from 
determining the radius during the ingress and egress. However, 
this effect is expected to be quite small. 



4. Parameter variation 

In order to determine the best fitting model parameters, a pa- 
rameter variation was performed to investigate the temperature 
and height dependence of the absorption lines. Furthermore the 
abundances of the three isotopologues found in the data was var- 



ied in order to infer isotopic ratios for '^C/'^C and "'O/'^O in 
Venus' atmosphere. 

Although a fixed temperature profile was used for the final 
fit, the profile was varied within the temperatures of the day 
and night side in order to quantify the response of the absorp- 
tion. Especially in the upper atmospheric layers above 100 km, 
large differences in temperature between the day and night side 
of Venus occur (see Fig.|4li, reaching about 150 K above altitudes 
of 170 km. Increasing the temperature from the night to the day 
side profile, the absorption of '^C '^Oa and '^C '^©2 increases by 
about 0.5% as can be seen in Fig.|5] The temperature variation 
has no significant effect on '^C "'O "^O. 




400 
Temperature [K] 



Fig. 4. Temperature profiles used for the calculation of transmis- 
sion spectra, obtained from lYung & DeMord (Il999h . Solid line: 
Day side. Dotted line: Night side. 



m 0.99 r: 




1597.6 1597.8 1598.0 1598.2 

Wavelength [nm] 

Fig. 5. Spectral response when varying atmospheric temperature 
profile. Solid line: Measured spectrum from the "Favorite-i-" re- 
gion. Dotted line: Spectrum calculated with the day side tem- 
peratures profile and Dash dotted line: Spectrum calculated with 
the night side temperature profile. Absorption line positions of 
the three isotopologues are indicated by different symbols. 



The tangential height of the beam was varied in steps of 1 km 
in the altitude range from 100 to 200 km. With an increase of 
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Table 1. Carbon and oxygen isotopic ratios 



Planet 



12/-/13, 



c/'-*c 



I6r,/lsi 



■o/'^^o 



Reference 



Earth 
Venus 
Venus 



89.01+0.38 
86±12 
110+30 



498.71+0.25 
500+80 

537+53 



(1) 
(2) 
(3) 



(n lYung & DeMora[T999l (2) lBezard etaPflgSTl (3) This work 



1 km in tangential height, the absorption of all lines decreases by 
about 1.3%. Since '-C "'O2 is the most abundant isotopologue 
and hence independent of the isotopic ratio, lines of this isotopo- 
logue are used to determine the best fitting tangential height of 
the measurements. Note that the determination of the tangential 
height is of course also dependent on the chosen beam width. 

Using initially Venus isotopic ratios of '^C/'^'C and '^0/'**0 
(see also Table [T]), isotopic ratios were varied by large amounts, 
i.e. in the range from half up to double isotopic ratios of Venus. 
The abundances of the CO2 isotopologues can be calculated ac- 
cording to the relationship from .Bezard et al., (.1987.) : 



12^/13^ 



'^COa/'^COa 



i^O/'^O = 2''c"'02/'-C"^0'^0. 



(1) 



It can be clearly seen that increasing the isotopic ratios, the abun- 
dance of '^C '*'Ot and '^c '^O "*0 decrease. 



5. Spectral fitting 

The fit of calculated and measured data was performed visu- 
ally since it was not possible to use a least-squares fit due to 
the poor spectral normalization and the remaining solar residu- 
als and noise features. Only isolated lines with a good baseline 
were used for the fitting which do not overlap with solar absorp- 
tion lines. 



5.1. "Favorite+" wavelength range 

In the "Favorite+" region (Fig.|3^), the '^C '^Oi absorption lines 
at 1 597.93 nm and at 1598.14nm (box symbols in the figure) 
can be used for the tangential height determination. The absorp- 
tion strength may however be disturbed by the higher baseline 
nearby. The best fit using these lines was achieved with a tangen- 
tial height of 1 16+1 km using a beam width of 15 km HWHM. 
Taking into account the uncertainty in the beam width (see sec- 
tion |2| the best fitting tangential height is 116+5 km. All other 
lines of this isotopologue overlap with other isotopologues or are 
at positions where solar absorption lines are located. Note that 
the inferred tangential height is close to the homopause located 
at about 130 km altitude, above which turbulent mixing of the 
atmosphere becomes unimportant and the atmosphere becomes 
stratified. Although this could influence the isotopic ratio over 
the altitude region probed, the atmospheric levels above the in- 
ferred tangential height do not contribute much to the absorption 
lines and this effect is neglected. Furthermore the data quality is 
too low as to resolve this effect. 

The absorption lines of '^C '^©2 at 1597.79nm and at 
1 598.08 nm (triangle symbols in the figure) are used for the 
determination of the carbon isotopic ratio, giving a best fit of 
'^C/'^C=11 0±30, which i s high er than the value for Venus in- 
ferred from iBezard et al.l (Il987l) . However taking into account 
the deviated basehne, the values are within the error margin 
comparable to literature values. Note that although the line at 



1 598.08 nm overlaps with '^C '^O '^O, the absorption is domi- 
nated by '^C '^02. 

The isolated lines of '^c isq "^O at 1 597.91 nm and 
1 598.27 nm (diamond symbols in the figure) could not be used 
for determining the isotopic ratio of '^0/'**0 due to the spike at 
about 1597. 89 nm. Nevertheless varying the '^0/"^0 ratio has 
only a marginal effect in the absorption. Therefore, the isotopic 
ratio found in the "Favorite" region (see hereafter) will be used 
here as well, with '^0/'^0=537±53. 

Figure |3^ shows the measured spectrum together with 
the best fitting model spectrum with a tangential height of 
116km, and isotopic ratios of '^c/i^c^llO and i60/''^0=537. 
Additionally, a scaled solar spectrum which has been extracted 
from one flatfielded image of the "Favorite+" region is plotted 
(dashed) to visualize the location of the solar absorption lines. 

5.2. "Favorite" wavelength range 

In the "Favorite" region (Fig. I3J)) the lines of '^C '^©2 are 
very weak and masked by absorption lines of the other iso- 
topologues and no isolated absorption lines can be found. In 
order to determine the best fitting tangential height, absorption 
lines of '■'C '^©2 are used, assuming that the isotopic ratio of 
'^C/'-'C is the same as in the "Favorite+" region. Absorption 
lines at 1 596.73 nm and 1 596.98 nm provide a best fitting tan- 
gential height of 120+2 km, when using a beam width of 24 km 
HWHM. Taking into account the uncertainty in the beam width, 
the estimated tangential height is 120+ 10 km. 

In order to determine the '^0/"^0 ratio, absorption lines of 
12c "'O '^O at 1 596.69 nm, 1 596.85 nm and 1 597.03 nm have 
been chosen. The first two lines are best fitted with '*0/''^0=575, 
whereas the other line is best fitted with '*O/"^O=500. Note 
that the line at 1596. 69 nm may be influenced by the strong 
baseline deviation shortward of this line and the overlap with 
a solar line. Therefore the estimated best fitting isotopic ratio is 
'^ 0/^^0=537+53. Thi s is also higher than the value determined 
bv lBezard et al.l (Il987h . but still in the range of their uncertainty. 

Figure |3j) shows the measured spectrum together with the 
best fitting model spectrum with a tangential height of 120 km, 
and isotopic ratios of '-C/'^C=1 10 and "'0/'**0=537. Like in the 
"Favorite+" region, a scaled solar spectrum is plotted (dashed) 
to visualize the location of the solar absorption lines. 



6. Summary and conclusion 

Using the Tenerife Infrared Polarimeter at the Vacuum Tower 
Telescope (VTT) in Tenerife, molecular absorption lines of the 
three most abundant CO2 isotopologues '-C "'O2, '''C '*02 and 
'^C '^O '^O have been detected within the atmosphere of Venus 
during its transit in June 2004. Although achieving limb spectra 
from the measured images was challenging, the spectra could 
be modeled using a line-by-line radiative transfer model, which 
even allowed to determine the isotopic ratios of '^c/'-^C and 
i6q^i8q^ which are slightly higher, but still in agreement with 
literature values. Note that the determination of isotopic ratios 
shall not serve as an additional measurement to what is already 
known about Venus, but rather as a proof of concept in order to 
show, what might be possible for exoplanets in the future. By 
furthermore varying the tangential height while leaving the tem- 
perature profile fixed, good agreement of modeled and measured 
spectra could be achieved with tangential heights in the range 
from (1 16±5)km to (120+10) km, depending on the wavelength 
range considered. 
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The observations resolved Venus' limb, showing an increas- 
ing Doppler-shift of the absorption lines with increasing altitude. 
This is probably caused by high altitude winds in Venus' at- 
mosphere, with a wind velocity gradient being consi stent with 
Venus Express measurements of lLellouch et alj (|2008|) . 

The results illustrate the feasibility of the transmission spec- 
troscopy technique to identify atmospheric constituents and even 
infer isotopic ratios in terrestrial exoplanets. Future planned 
space missions, which are designed to measure the atmospheres 
of terrestrial exoplanets, will utilize similar techniques to exam- 
ine their atmospheres in detail. 

7. Implications for extrasolar planets 

Current models for exoplanets mostly rely on knowledge derived 
from the solar system planets. It is straightforward (at least con- 
ceptually) to build a spectrum from a given planet type. The in- 
verse problem, however i.e. to infer the characteristics of a planet 
from a spectrum is much more difficult due to its ill-posed na- 
ture. 

Transiting exoplanets cannot be resolved in front of their 
central star as was the case for the Venus transit presented here. 
Thus the method to obtain atmospheric transmission spectra of 
an exoplanet is different than it was performed in this work. For 
transits of exoplanets, where the planet and its central star are 
convolved within a few pixels, the total flux of the system has to 
be compared during and outside the transit. Since the apparent 
radius of the planet and thus the transit depth of the transiting 
planet is dependent on the wavelength due to the absorbing at- 
mosphere, a transmission spectrum can be obtained by observing 
the transit at diff'erent wavelengths. 

Current and near-future space telescopes provide only low 
resolution spectroscopy (Spitzer: R<600, James Webb Space 
Telescope: R<3000), thus only molecular absorption bands can 
be investigated rather than single absorption lines. Although 
ground-based telescopes already feature high-resolution spec- 
trographs of R ~ 100 000 like the CRyo genic high-resolut ion 
Infrared Echelle Spectrograph (CRIRES, JKaeufl et al.ll2004 at 
the Very Large Telescope (VLT), the SNR is still too low to 
detect lines of different isotopologues as was performed here, 
since atmospheric absorption lines of exoplanets are contrasted 
against the overall stellar flux. Atmospheric features of exoplan- 
ets are in the order of magnitude of around 10"^, which is much 
smaller than the 1 to 2% absorption detected for Venus here. This 
makes the investigation of the atmosphere much more challeng- 
ing. Furthermore, ground-based observations have to remove tel- 
luric absorption lines in order to access the atmospheric absorp- 
tion of an extrasolar planet. 

Nevertheless, the investigation of an terrestrial exoplanet 
is alr eady possible, as was shown for GJ 1214b (iBean et alJ 
I2OIOI) . The data can however be interpreted in several ways 
and the actual compositi on of the atmosphere is still sub- 
ject to discussion (see e.g. ICroll et al]l201 it iDesert et al]l201 it 



The thermal emission of about 30 exoplanets has already been 
detected, which allow the determination of the temperature 
structure, energy redistribution and day-night contrasts (se e e.g. 



iMiUer-Ricci Kempton et al.ll201 ll) 



Venus composition and p, T profiles are well known from 
satellite and ground-based observations, which is of course 
not the case for exoplanets. Temperature profiles of exoplanets 
might be inferred by emission spectroscopy during secondary 
eclipse observations. For example, CO2 absorption bands seen 
in emission would indicate a temperature inversion in the strato- 
sphere since CO2 can be assumed to be well mixed through- 
out the atmosphere, as is the case for Venus, Mars and Earth. 
Also the surface temperature might be inferred in spectral win- 
dows where the atmosphere is transparent down to the ground. 



iDeming et al.]l2005l: [Harrington et alJl200alKnutsonll2007h . 

The determination of exoplanet isotopic ratios (as shown in 
this paper for Venus) would give insight into the evolution and 
diversity of exoplanet systems. However this is beyond the abil- 
ity of current and near future observational efforts, because both 
a significant SNR and a high spectral resolution are required to 
detect and separate the lines: 

The Venus transit was observed with a spectral resolution of 
R =200 000 in a wavelength range, where the three most abun- 
dant isotopologues are detectable. In the wavelength range ob- 
served, a resolution of more than about 25 000 is required in or- 
der to separate lines of '^C '^Oi and '^C '^©2 to determine the 
'^C/'^C ratio. For the detection of '-C '*0 '^O absorption fines 
for the determination of the "'0/"*0 ratio, a much higher spec- 
tral resolution of about R =50 000 is required, to separate them 
from the '^C ^^02 and '^C '^02 absorption lines. However, note 
that other molecules might feature stronger absorption lines that 
are well separated. 

iRauer et aP (1201 ih have shown for the transmission spec- 
trum of an Earth-like planet around a solar like G star at 
a distance of lOpc that even with a low spectral resolution 
of R =2000 the photon-limited SNR is about 0.054 for the 
4200nm CO2 band, when using e.g. the James Webb Space 
Telescope (JWST) having an aperture of 6.5 m. At R =10 000 the 
corresponding SNR would be 0.022 and at R =50000 it would 
be 0.009, hence hundreds of transits need to be co-added in order 
to achieve a reasonable SNR of more than three. Since the SNR 
is directly proportional to the telescope aperture, a larger tele- 
scope aperture would increase the SNR. However, only with next 
generation telescopes having much bigger apertures of about 
40 m (e.g. the European Extremely Large Telescope, E-ELT), 
SNRs close to unity are achievable for terrestrial exoplanets. 
Observations of close-in hot-Jupiter planets will also signifi- 
cantly increase the SNR. Thus the determination of isotopic ra- 
tios of an exoplanet atmosphere is expected in the near future. 

Despite the difficulties for high-resolution spectroscopy of 
extrasolar planets, the observations of Venus shown here illus- 
trate the potential of such analysis. Once sufficiently sensitive 
instrumentation becomes available, the search for isotopic lines 
in exoplanet atmospheres can provide severe constraints on the 
atmospheric evolution of these distant planets. 
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